Few-body spin couplings and their implications for universal quantum computation 
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Electron spins in semiconductor quantum dots are promising candidates for the experimental 
realization of solid-state qubits. We analyze the dynamics of a system of three qubits arranged 
in a linear geometry and a system of four qubits arranged in a square geometry. Calculations are 
performed for several quantum dot confining potentials. In the three-qubit case, three-body effects 
are identified that have an important quantitative influence upon quantum computation. In the 
four-qubit case, the full Hamiltonian is found to include both three-body and four-body interactions 
that significantly influence the dynamics in physically relevant parameter regimes. We consider the 
implications of these results for the encoded universality paradigm applied to the four-electron qubit 
code; in particular, we consider what is required to circumvent the four-body effects in an encoded 
system (four spins per encoded qubit) by the appropriate tuning of experimental parameters. 
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I. INTRODUCTION 

Electron spins in semiconductor quantum dots are a 
leading candidate for the physical realization of qubits 
in a quantum computer Although any quantum al- 
gorithm can be implemented using single-qubit and two- 
qubit gates 0, many such algorithms realize substantial 
increases in efficiency by exploiting simultaneous interac- 
tions among three or more qubits [E EL IE IE II IE IE EE 
EE O EE El EE 111 • In order to employ such simulta- 
neous interactions, it is essential to understand in detail 
the many-body dynamics of the system of coupled qubits. 
More generally, since a practical quantum computer may 
need to contain as many as 10 6 qubits it is essential 
to characterize the effect of many-body interactions on 
the system's overall energy landscape. 

In past work 0,0, we used a model confining poten- 
tial of superposed parabolic minima to demonstrate that 
three-body effects significantly influence the Hamiltonian 
of three electrons confined to three quantum dots at the 
vertices of an equilateral triangle and that four-body ef- 
fects are significant for four electrons confined to a tetra- 
hedral arrangement of four dots. Here we extend these 
results in two ways. First, we analyze three quantum dots 
in a linear geometry |19| and four dots in a square geome- 
try |20j | since these geometries are more likely to occur in 
a real quantum computer apparatus. Second, by employ- 
ing a Gaussian shape for the confining potential of each 
well [2EI , we explore the sensitivity of the many-body ef- 
fects to the form of the confining potential. In both cases, 
a non-perturbative calculation finds that many-body ef- 
fects contribute appreciably to the Hamiltonian. We note 
that Scarola et al. |22j, |2jj have demonstrated that the 
application of a magnetic field allows chiral terms to arise 
in the spin Hamiltonian, which modifies this Hamiltonian 
in another important manner as compared to the naive 
Heisenberg form. 

To date, discussions of quantum dot quantum compu- 



tation have nearly always assumed pairwise Heisenberg 
interactions. In view of the above result, this implies that 
computational errors may occur in the context of quan- 
tum computers using electron spin qubits in quantum 
dots, unless one always simultaneously couples only dis- 
joint pairs of dots. There are at least four circumstances 
where this may be undesirable or even infeasible. One is 
fault tolerant quantum error correction, where simulta- 
neous operations on several coupled dots have been asso- 
ciated with better error thresholds. A second is adiabatic 
quantum computation |Tlj . in which the final Hamilto- 
nian may include the simultaneous interactions that we 
discuss here. We will not analyze these possibilities here, 
although we believe that the methods we discuss below 
are relevant to them. 

We will focus on two other contexts, that of "encoded 
universality" (EU) [E IE II IE H and that of compu- 
tation on decoherence-free subspaces (DFSs) [E IE EH 
and supercoherent qubits [l(|. In these cases, the goal is 
to perform universal quantum computation using (EU:) 
only the most easily controllable interaction, or (DFS, 
supercoherence:) using only interactions that preserve 
the code subspace, since that subspace offers protection 
against certain types of decoherence. (Strong and fast ex- 
change interaction pulses can further be used to suppress 
decoherence [2^] and to eliminate decoherence-induced 
leakage prjj.) 

We will refer to these cases collectively as "encoded 
quantum computation." It turns out that universal 
quantum computation using only the Heisenberg ex- 
change interaction is an extremely attractive possibility 
in encoded models, and we will consider it in detail be- 
low. After establishing that four-body interaction terms 
can arise in a Heisenberg exchange Hamiltonian, we in- 
vestigate the question of neutralizing their effect by using 
encoded qubits [E IE IE IE EE SEE EE El EE S3 . By 
generalizing the work of Bacon |30j , who showed that uni- 
versal quantum computation was possible using encoded 
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gates with two-body coupling Hamiltonians (i.e., assum- 
ing that the Heisenberg Hamiltonian was applicable even 
when coupling three or more dots at a time), we enumer- 
ate tuning conditions on experimental parameters that 
are needed for the four-body effects to cancel out. An al- 
ternative is to design these encoded gates while allowing 
only pairs of electrons to couple at any given time. This 
is indeed possible, as shown in Ref. [3l]], for the price 
of significantly longer pulse sequences per given encoded 
gate. Nevertheless, in view of the findings reported here 
and in Refs. [22, |23| , this price may be worth paying. 



II. THREE-ELECTRON CASE 

A system of three electrons within a confining scalar 
potential V(r) obeys the Hamiltonian 



eight (unnormalized) basis states is 
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in the absence of spin-orbit coupling and external mag- 
netic fields. Although Ref. [TtJ succeeded in demon- 
strating significant three-body and four-body effects in 
systems containing three or more electrons, a confining 
potential with quadratic minima has certain other char- 
acteristics which are unlikely to describe an experimental 
arrangement; for example, it diverges at large distances 
from the origin, and the single adjustable parameter to 
forces us to specify very narrow minima whenever we 
want a high barrier between them. We therefore begin 
with the Gaussian form 



V(v) 
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which has two tunable parameters. The three fixed 
points are collinear and separated by a distance 21: 
A = (-21, 0, 0), B = (0, 0, 0), and C = (21, 0, 0). 

We assume a Heitler-London approximation [32l| . 
wherein excited orbital states and states with double oc- 
cupation of any single dot are neglected (see Ref. [23j for 
a recent discussion of the validity of this approximation 
in the context of electron spin qubits). The system's only 
degrees of freedom are therefore the spins of the confined 
electrons, leading to a total of 2 3 — 8 "computational" 
basis states 

\^(s A ,s B , Sc )) =Y,S P P[\A)\B)\C)\ sa) \sb) \sc)]- 
p 

(4) 

In the above, \{A}) are the three localized orbital ground 
states; denote the corresponding spin states; P is 

the set of all permutations of {^4, B, C}; and 5p is 1 (-1) 
for even (odd) permutations. For instance, one of the 
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To characterize the localized orbital state \{A}) for each 
dot, we expand @ to quadratic order and solve the 
Schrodinger equation as though the other potential wells 
were absent: 



fmuj \ 3 / 4 ( muj . 
(r) = (r\A) = (_ ) exp(- — |r- A|" 



(5) 



Unless a is small compared to l~ 2 , of course, this is a 
much coarser approximation than it would be for purely 
quadratic minima, so we refine it by centering 4>a (r) and 
4>c( r ) at the points which minimize and (C\h\C). 

Because these orbitals overlap at least slightly for any 
finite u! a , the states Q are not orthogonal. 

We now define H sp i n to be the matrix representation 
of H in the basis l@J), and expand it in terms of tensor 
products of Pauli matrices: 



This expansion is always possible, since the set of n- 
fold tensor products of Pauli matrices constitutes a com- 
plete orthonormal basis for the linear vector space of 
all 2™ x 2™ matrices. Because we have written the ba- 
sis (0J in the form \sa) \sb) \sc), these Pauli matri- 
ces can be associated with spin operators on each of 
the three quantum dots. For example, we can write 
g\ ® 0-3 ® oo = 2S a ,x® c 2Sb,z®I = ^Sa, x Sb,z, where the 
notation Sw,i means the Pauli operator Ui applied to the 
electron in the quantum dot at W , and where / is the 
2x2 identity matrix. (We exclude fi from the definition 
of the matrices af, thus, the c^fc have the dimensions of 
energy.) In the case of an arbitrary 8x8 matrix, 64 com- 
plex numbers would be required to specify our , but 
the operator clearly has certain properties which con- 
strain the values of the coefficients, such as Hermiticity, 
reflection symmetry, rotation symmetry, inversion sym- 
metry, and invariance under permutation of the electrons' 
labels. Once these symmetries have been accounted for, 
the Cijk may be characterized by just three real quanti- 
ties: 

H sp in - K +K 2 [AB}(S a -Sb + Sb-Sc)+K 2 [AC]Sa-Sc, 

(6) 

where Sw ■ = Sw x Sy tX + Sw^ySy.y ~\~ Sw z Sv,zi and 
1^2 [A?] i s the pairwise coupling coefficient between the 
spins of the electrons in dots i and j. Here and else- 
where, we use symmetry considerations to reduce the 
number of coupling coefficients in our equations; in this 
case, the reflection symmetry of (0 through the x-z plane 
implies that K 2 [AB] = K 2 [BC]. Physically, the constant 
K 2 [AB] quantifies the coupling between adjacent spins, 
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while K 2 [AC] describes the coupling between the spins 
at opposite ends of the row. 

Defining St = ^>a + S_b + Sc, one finds that 



H. 



spin 



where 
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2L X 

2Li + 2L' v 



where the indices on the left-hand side denote the val- 
ues of St, St, z , and |S^ + Sc| respectively. Although 
the states | ^{sa, sb,sc)) are not orthonormal, the eight 
states (0 are orthogonal, and they are also eigenvectors 

(7) of the 8x8 matrix J7J), which means that H sp i n has been 
diagonalized. To obtain the parameters {£n, L\, L^}, we 
will choose three eigenstates with different good quantum 
numbers, and observe that their energies can be evalu- 

(8) ated either by matrix algebra or by integrating micro- 
scopically over the axes and the spins to compute the 
expectation value of (|TJ: 



The expansion JJJ reveals that any simultaneous eigen- 
state of (Sa + Sc) 2 and S T is also an eigenstate of H sp i n . 
We can construct such simultaneous eigenstates by us- 
ing the Clebsch-Gordan table twice, first to combine the 
spin of the electron in dot A with the spin of the electron 
in dot C, and then to combine that spin-1 (or spin-0) 
system with the spin of the electron in dot B: 
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Inserting (JJJ into the left-hand side, for three distinct 
combinations of the good quantum numbers {(S^ + 
S c f,S T }, yields 
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Lo + —Li + 2L\ (11) 
L + jLi + 2L' X 



^0 + ^1, 



while the corresponding wave functions © turn the 
right-hand side into 



(tt(ttt)|g|tt(ttt)) 

(*(TTT)|*(TTT)> 

<^(Tn)|g|^(tn)) + 2(w(nt)|g|^(nt)) - 4(^(rn)|g|^(nt)) + wpto)) 

(*(TTI)I*(TTI)> + 2(*(tit)l*(TiT)> - 4(*(TTI)I*(TIT)> + (*(m)|*(m)> 

(^an)|g|^(Tn)) - (garrjjgjgurr)) 
(*(m)i*(m)> - <*(m)i*an)) ' 



The evaluation of these matrix elements and overlap in- 
tegrals is a tedious, but straightforward procedure given 
the microscopic forms of H and tp(r) in QJ, Q, and JSJ. 
Combining Eqs. (JTIJ, and (|T2l . we thus compute 

K , K 2 [AB] , and A 2 [AC] in terms of oj and the dimen- 
sionless system parameters 



2,2 



muj„l 



Xb 



„2 



nlhuj 
2Vp 

huj 



(13) 
(14) 
(15) 



Physically, the quantity is the ratio of the height of 



the potential barrier between wells to the energy of the 
orbital ground state JSJ, while x c is the ratio of the equi- 
librium Coulomb repulsion potential to the energy of the 
orbital ground state, and x v is the ratio of the individual 
well depth Vq to the ground state energy. 

Here and in the following section, we have estimated 
experimentally relevant values of Xb and done in 

Ref. [l| . We assume that the width of the function J5j , 
which is 2y/ h/mcu , must be roughly equal to the sepa- 
ration between adjacent dots 21; using l|13f) . we conclude 
that Xb 1. For GaAs heterostructurc single dots, k ss 
13, m* rs 0.067 m c , and Hu « 3 meV, which according 
to l|14f> means that x c w 1.5. 

A potential of the form © is most suitable for quan- 
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turn computation when al 2 is close to 1; if the inverted 
Gaussian decays too quickly in space, the spin coupling 
in the system becomes negligible, and if it decays too 
slowly, the local minima in V tend to coalesce at the 
center. Using ±Tiuj ~ 1 meV, Vo ~ 3 meV 0, and 
our prior estimate of x c ~ 1.5, we obtain the relation 
xi, Ri x v ~ 3, by applying ljT3|) . ((Tl*|l . and i(T5|) . Noting 
that the parameter x c has very little influence on any 
of the coupling constants over physically realistic ranges 
of xi, and x v (and in any event depends on quantities, 
such as k, which are difficult to tune experimentally), we 
henceforth set x c — 1.5. 

Fig. n shows the energy shift K as a function of 
the system parameters {xb,x v }. As one might expect, 
this spin-independent quantity increases with increasing 
x v and decreasing Xf, (whenever lo d decreases, there is 
greater orbital overlap and thus more Coulomb repul- 
sion, irrespective of spin state) . The coupling constants 
K 2 [AB] and K 2 [AC] are plotted in Fig. [land Fig. re- 
spectively. We notice that they differ (which rules out 
the simple Heisenberg form H spin = JJ2i<j(^i ' ^j)), 
and that K2[AC] is only about an order of magnitude 
smaller than K 2 [AB], as we have confirmed by study- 
ing K 2 [AB](xb, x v ) and K 2 [AC](xb, x v ) on a logarith- 
mic scale. In the context of quantum computation, this 
demonstrates that a nearest-neighbor approximation for 
the coupling between dots is insufficient (see also Ref. 
[23| , where a similar conclusion was reported using a low- 
energy Hubbard model with one electron per site). 



III. FOUR-ELECTRON CASE 

For the case of four quantum dots arranged in a square 
of side 21, our formalism is more complex in detail but 
identical in structure. We therefore describe the compu- 
tation only in outline. 

The confining potential in the coordinate Hamiltonian 



Expanding H in terms of products of Pauli matrices, 



i=l 

now becomes 
V(r)-- 
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-a|r-B| 2 
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where A = (0,21,0), B = (2Z,2Z,0), C = (2Z,0,0), and 
D = (0,0,0). Our computational basis consists of 16 
fully antisymmetrized vectors of the form 

\y( SA ,SB,sc,s D )) = Y / ^P p U)\B)\C)\D) 
p 

®\sa)\8b)\8 C )\8d)]. (17) 

The form of <p(r) remains the same; to maintain the re- 
quired geometrical symmetries, we now shift all four lo- 
calized orbital wave functions an equal distance toward 
the point (1,1,0). 
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we discover by applying the symmetries of (|16fl that four- 
body terms now appear with nonzero coupling coeffi- 
cients: 

#spin = K Q + K 2 [AB](S A -Sb + Sb-Sc 
+Sc ■ Sd + Sd ■ S A ) 
+ K 2 [AC](S A ■ S c + S B ■ S D ) 
+ K 4 [ABCD][{S A -S B ){Sc-S D ) 

+ (S B -S C )(S D -S A )} 
+ Ki[ACBD]{S A ■ S C )(S B ■ S D ), (18) 

where K^[ijkP\ is the four-body coupling coefficient 
among the spins of the electrons in dots i, j, k, and 
I. Physically, the constant K 2 [AB) describes the pair- 
wise coupling between adjacent spins, while K 2 [AC] de- 
scribes the pairwise coupling between non-adjacent spins, 
Ki[ABCD\ describes four-body interactions concentrat- 
ing on pairs of adjacent spins, and K^ACBD] describes 
four-body interactions concentrating on pairs of non- 
adjacent spins. We define St — S A + Sb + Sc + Sd, 
which leads us to 



H, 



spin 



where 



Lq + L X S| + L[[(S A + Sc) 2 + (S B + S D ) 2 } 
+L 2 (S 2 T ) 2 + L' 2 (S A + S C ) 2 (S B + S D ) 2 , (19) 
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K A [ABCD] 
KAACBD] 



45 

L + 3Li + 3L[ + —L 2 - 
2Li + 24L 2 

2ii + 2L[ + 2AL 2 + 3L' 2 
8L 2 

8L 2 + 4L 2 . 



-L' 2 (20) 



Applying the Clebsch-Gordan table three times creates 
sixteen simultaneous eigenstates of (S A + Sc) 2 , (Sb + 
Sd) 2 , and S T . Inserting five of these states with dif- 
ferent quantum numbers into l|10|) yields five equations 
for the five unknowns {Lq, L%, L[, L 2 , L 2 } in terms of the 
eigenstate energies. As before, these energies may be ex- 
pressed in closed form as functions of Xb, x c , and x v by 
integrating the right-hand side of l|10(l explicitly. 

The energy shift K for the square case is plot- 
ted in Fig. 01 as before, this constant is largest for 
strongly Coulomb-coupled dots separated by low po- 
tential barriers. Figs. 03 H3 an d El depict the 
coupling coefficients K 2 [AB], K 2 [AC], K 4 [ABCD], and 
Ki[ACBD] respectively. The departure from the pair- 
wise Heisenberg picture is even more pronounced here: 
we see that for physically relevant values of the param- 
eters {xb,x v }, the four-body coefficient K^ACBD] is 
of the same order of magnitude as the two-body coef- 
ficient K 2 [AC], while K 4 [ABCD]/K 2 [AB] - 0.1, as is 
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confirmed by plotting K 2 [AB](xb,x v ), K 2 [AC](xb, x v ), 
K^[ABCD](xb,x v ), and Ki[ACBD]{xb, x v ) on a loga- 
rithmic scale.. Typically, K^ACBD] is opposite in sign 
to K 2 [AC] , leading to a particularly important competi- 
tion between the two-body and four-body interactions. 

In order to confirm that the qualitative similarities be- 
tween our final results and those of Ref. 01 were not ar- 
tifacts of having made two broad changes to V(r) rather 
than one, we also analyzed both the N = 3 and N = 4 
dot geometries using a confining potential of superposed 
quadratic minima. The variation of the coupling coeffi- 
cients, within experimentally relevant ranges of Xb and 
x c (analogous to Figs. ^ through [SJ, strongly resembled 
that for the Gaussian potential in all cases. 



i, j be denoted as 
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Then a single encoded DFS qubit is formed by the two 
singlets of four spins, i.e., the two states with zero to- 
tal spin St = | Sa + + Sc + S^ | . These states are 
formed by combining two singlets of two pairs of spins 
(|0l)), or triplets of two pairs of spins, with appropriate 
Clebsch-Gordan coefficients (|1l)): 



IV. COMPUTING IN THE PRESENCE OF 
FOUR-BODY INTERACTIONS USING 
ENCODED QUBITS 

We have shown that coupling three dots simultane- 
ously quantitatively modifies the value of the exchange 
constant, and that coupling four dots simultaneously 
switches on a four-body interaction term of the form 
Ki[ABCD](S A ■ S B ){S C • Su) and its permutations. 
This conclusion appears to be robust under changes in 
dot geometry and in the confining potential. A natu- 
ral question is whether there exist methods to cancel 
the four-body correction. The issue is particularly ur- 
gent when one considers encoded quantum computation 
(EQC). In many known constructions of universal gates 

for eqc : i i 1 i I M M El ES El M ES & 

there arises the need to simultaneously couple several 
spins. One of the most popular codes, described in de- 
tail below, uses four spins per encoded, or logical qubit 
II IS II IS IIS EE M El El El For this code, uni- 
versal computation requires that four spins be coupled 
at the same time using pairwise Heisenberg interactions. 
Hence a priori it appears that EQC using the four-qubit 
code suffers from a fundamental flaw. We now explore 
whether the four-qubit code may be implemented in such 
a way that each four-body coupling is either cancelled or 
reduced to an overall phase. Our findings highlight prob- 
lems that the four-body terms present in the context of 
EQC, and also provide an interesting perspective on how 
the four-body terms may need to be dealt with in general. 



The Code 



Let us descibe the four-spin DFS code, first proposed 
in Ref. |27| in the context of providing immunity against 
collective decoherence processes (see Ref. [2i| for a re- 
view). Let the singlet and triplet states of two electrons 



|0 L ) = \s)ab ® \s)cd 

1 



^(i*(tm)> + i*(mr)) 

-|*(TUT)>-|*UtU)» 



(21) 
(22) 



= 



1 



-j=(\t-) AB ®\t+)cD 
+\t+)AB®\t-)cD) 



\to)AB ® |*0>CZ» 



= _(2|*(TTU)) + 2|*(UTT)> - l*(tUT)> 

-i*um)> - i*(uu)> - i*uut)». (23) 

As shown in Refs. [S |(J , the Heisenberg interaction Si • Sj 
can be used all by itself to implement universal quantum 
computation on this type of code. The Heisenberg in- 
teraction is closely related to the exchange operator By , 
defined as 



1 



EL 



1 

1 



(24) 



via Eij — | (4Sj ■ Sj + 1) . The difference in their action as 
gates is only a phase, so that we will use and S, ■ Sj 
interchangeably from now on and write Eij ~ Sj • Sj . The 
have a simple action on the electronic spin up/down 
states, as seen from the matrix representation (|24|) : the 
states 1 00) and |11) are invariant, whereas |01) and 1 10} 
are exchanged. Using this, it is simple to show that, 
in the {|0l), |10} basis, the exchange operators can be 
written as 6, 28] 



Eab — Ecd — 



-1 
1 



= -Z 



(25) 



Eac — Ebd — + —Z 

Ead — Ebc = —^-X + -Z, 

where X, Z are the encoded Pauli matrices cr x ,a z , 

i.e., the Pauli matrices acting on the |0l),|1l) states. 
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It follows from the Euler angle formula, e ~ zujn cr = 
e -i@<rz e -i0°-x e -ia«Tz ( a ro t a tion by angle lj about the axis 
n, given in terms of three successive rotations about the 
z and x axes), that one can perform all single encoded- 
qubit operations on the DFS states, simply by switching 
the exchange interaction on and off. Note that the Eu- 
ler angle formula is satisfied by any pair of non-parallel 
axes, but orthogonal axes may be more convenient. One 
can obtain an encoded a x operation by switching on two 
interactions simultaneously for the appropriate time in- 
tervals: 

X = -2 (e ac + \e A i^J /\/3 = (E AC - E AD ) /Vs. 

Use of the Euler angle formula requires a Hamiltonian 
which is a sum of exchange terms with controllable coef- 
ficients Jij{t): 

H s = ^2 Jij(t)Eij. 

This is achievable, e.g., by using local magnetic fields 
Q ED, EJ IHiby ferroelectric gates ^3, or by opti- 
cal rectification |3fij . It is important to emphasize that 
the last two methods |H |3g do not require magnetic 
field control, hence overcome at least in part the prob- 
lems with EQC raised in Refs. jHH^]. This is an impor- 
tant advantage with regards to EQC, which renders these 
electrical-only type control methods distinctly preferable 
to those using magnetic fields. However, residual mag- 
netic fields, e.g., due to nuclear spin impurities, do remain 
a problem, especially in the group III-V semiconductors, 
such as GaAs [37$. In silicon-based architectures this 
problem can be minimized by isotopic purification |3S| ]. 

B. Effect of the Four-Body Terms on a Single 
Encoded Qubit 

Let us now consider how the four-body terms act on 
the DFS code. Using the results above, we find that 

{S A • S B )(S C • S D ) ~ EabEcd = {-Zf = I, 



where / is the identity operator. Also, 



E A cEbd — ^~2~^ 2 y 

= -(3I + I+ V3(XZ + ZX)^ = I, 

and similarly E A jjEbc = I- Thus all fourth-order terms 
(Si ■ Sj)(Sfc - Si) oc I as long we restrict their action to 
the subspace encoding one qubit. This implies that the 
encoding into the 4-qubit DFS is immune to the fourth- 
order terms. In other words, when this encoding is used, 
the problem of the computational errors induced by the 
undesired fourth-order terms simply disappears, as long 
as we restrict our attention to a single encoded qubit. 
C. Two Encoded Qubits 

We must also be able to couple encoded qubits via 
a non-trivial gate such as controllcd-phasc: CP = 
diag(— 1, 1, 1, 1). This is one way to satisfy the require- 
ments for universal quantum computation |3£} , though it 
is also possible to complete the set of single-qubit gates 
by measurements |40j . Two encoded qubits of the form 
(12211 . (|23fl occupy a four-dimensional subspace of the zero 
total spin subspace of 8 spins. The zero total spin sub- 
space is 14-dimensional. A very useful graphical way of 
seeing this, introduced in Ref. Q but also known as a 
Bratteli diagram, is depicted in Fig. [5] 

As more spins are added (horizontal axis), there are 
more possibilities for constructing a state with given total 
spin (vertical axis). In the case of four spins there are two 
paths leading from the origin to St = 0; these correspond 
exactly to the |0l) and code states. For eight spins 
there are 14 such paths. Only four of these correspond to 
the four basis states |0l0l), |0l1l), It is 

convenient to label paths according to the intermediate 
total spin: the state \S\, S2, S3, S4, S5, Sq, S7, Sg), where 
Sfe is the total spin of k spin-1/2 particles, uniquely cor- 
responds to a path in Fig. (we omit the origin in this 
notation), and the Sfc form a complete set of commuting 
observables 0- E.g., 



|0 L L ) = |l/2,0,l/2,0,l/2,0,l/2,0) 
|0 L 1 L ) = |l/2,0,l/2,0,l/2,l,l/2,0) 

|U0i) = 11/2,1,1/2,0,1/2,0,1/2,0) 

|1 L 1 L ) = |l/2,l,l/2,0,l/2,l,l/2,0) 

I — 



/ \ / \ / \ 

/ \ 

/ \ / \ / \ 

/ \ / \ 

/ \ / V 



On the right we have indicated the path in Fig. cor- responding to each state. The other 10 states with 
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zero total spin can be similarly described. Thus the set E.g., in this basis the operator Ede has the 
of 14 states {\S\, S%, S3, S4, £5, 5*6, £7, 0)} forms a basis representation |4lj 
for the subspace of zero total spin of 8 spin- 1/2 par- 
ticles. Henceforth we will find it convenient to repre- 
sent exchange operators in this basis. We will order 
the 14 basis states as follows: first the four code states 
\0l0l), \0lU), |U0 l ), \1 l 1 l ) as above, then 

|l/2,0, 1/2, 1, 1/2,0, 1/2,0), |l/2, 1, 1/2, 1, 1/2,0, 1/2,0), 
|l/2,0, 1/2, 1, 1/2, 1, 1/2,0), [1/2, 1, 1/2, 1, 1/2, 1, 1/2,0), 
|l/2,0, 1/2, 1, 3/2, 1, 1/2,0), |l/2, 1, 3/2, 1, 1/2, 0, 1/2,0), 
|l/2, 1, 3/2, 1, 3/2, 1, 1/2, 0), [1/2, 1, 3/2, 2, 3/2, 1, 1/2, 0), 
[1/2, 1, 1/2, 1, 3/2, 1, 1/2,0), [1/2, 1, 3/2, 1, 1/2, 1, 1/2,0). 



2 



E 



DE 



V3 
2 



2 



2 



VI 

2 



2 



2 



2 



4 

V^L5 
4 



_15 

4 

_L 

4 



1/ 



Recall that the first four rows refer to the code space. It 
is then clear that Ede mixes the code space with four 
of the remaining ten states that have zero total spin. 
This is a general feature of all exchange operators acting 
on two code blocks simultaneously. For this reason it is 
impossible to couple two code blocks in one step, while 
preserving the code space. 



to the appearance of three- and four-body corrections. 
For ease of visualization, we will assume that the two 
blocks are squares of side 21 and that dots D and E are 
separated by a distance 21, although nearly all of the fol- 
lowing calculations are independent of the exact spatial 
relationship between the blocks. 

The gates are (adapting the notation of |3fJ, App. E[) 



D. Enacting an Encoded Controlled-Phase Gate 

For the 4-qubit code above, procedures implementing 
a CP gate were first given in Refs. |5|,|6|. Recently, Ba- 
con [30j, App. E] found a simplified scheme which is a 
useful starting point for our purposes. Bacon's imple- 
mentation of a CP gate between two pairs of 4-qubit 
blocks (qubits A-D and qubits E-H) involves a sequence 
of 14 elementary gates, each of which requires that sev- 
eral simultaneous exchange interactions be switched on 
and off. We will here take the approach of utilizing Ba- 
con's construction, while making some modifications due 



u 3 
u 6 



exp 

exp 
exp 
exp 



' • / 1 D \ 

— j= ^~~3Eef — - {Efg + Efh + Eqh) 
in ( _„ 2 



- \ SEcd — - (Eab + Eac + Ebc) 
4V2 V 

i| (e fg + \e gh 



(u A u B u\ui 



2 

IJI > 
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where 



Ua — exp 
Ub = exp 



■cos" 1 {-1/3)E DE 



ITT ° 

A=Kj 



In terms of these gates, the controlled-phase gate can be 
written as 

CP = ul(U^ul)U^U 6 U 6 (PsU a )Ui, 

where U3U2 can be executed in one step since by in- 
spection the two gates operate on the two blocks sepa- 
rately (and identically). This gate sequence operates in 
the entire 14-dimensional subspace of St = states of 8 
spins: the code space is left after application of U\, but 
is returned to at the end of the sequence, when u\ is 
applied. Hence our single-qubit considerations above do 
not apply: even if a four-body interaction acts as the 
identity operator on a single encoded qubit, it may act 
non-trivially in the larger St = space. We must there- 
fore carefully analyze the action of this gate sequence in 
light of the three- and four-body corrections. 



The key point in Bacon's construction of the gate se- 
quence is to ensure that each gate acts "classically," i.e., 
it only couples a given St = basis state to another, 
without creating superpositions of such basis states (that 
the gates above act in this manner is not at all simple 
to see directly, but is the reason for the particular choice 
of angles in the gates). Here we will show that, in order 
to still satisfy this key point, it is necessary to tune the 
four-body exchange coupling constants. Thus, to enact a 
CP gate in the presence of four-body interactions, there 
needs to be sufficient flexibility in tuning the four-body 
coupling. We note that there are other ways to obtain 
a CP gate |3l): our point here is mostly to explore the 
implications of the four-body terms in a context of some 
general interest. Let us now consider each of the gates in 
detail, in increasing order of complexity. 



1. The Ua Gate 

Ua only involves a single exchange interaction, and so 
is unmodified in the presence of the three- and four-body 
corrections: 



U' A 



Ua = exp 
/ a 



cos- 1 (~l/3)E DE 



1 

iV2 



1 

i73 



a 



1 

1V2 



1 

iV2 



1 

iV2 



1 

iV2 



1 

iV2 



1 

iV2 



ft 



(3 



1 1 1 J_ /5 

73 2iV6 2» V 2 

j_ IT _i 1 

2-i V 2 V3 2iV5 



(3 



where a = and j3 = e -V(2co S I (-3))_ 



2. The U B Gate 



three-body correction to the exchange constants. In ad- 
dition, a coupling between dot F and dot H will arise, 
which forces our modified U5 gate to have the form 



U5 involves dots F, G, and H, with dots F and H 
simultaneously coupled to G, and hence experiences a 
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= exp 



( p $ 

$ P* 



p $ 



t(2J 5 +3) 
2V3 



iT(2J 5 +3) 
e 2 v^ 



i-ir(2J 5 +3) 

e 



t(2J 5 +3) 
2\/3 



where 



3. The U B Gate 



Bacon's Ub gate is 



A 

V 



v /|(4) 2 -24 + i 



cos 



-A 



V3A 



-(1-4) sin 



-A 

2 . 



-A 



It is seen that this gate operates "classically," and is non- 
diagonal, only when J' 5 is either or chosen such that A 
is an even integer. In both cases, we recover Bacon's 
functional form: 



i 



U 5 



- V3ir 

e l 2 



e l 2 



„■ V3Tr 

e 2 



e* 2 



exp 



51 Ei 3 



K A=i<j 

= diag(l, 1,1, 1,-1, -1,-1, 

-1,-1,-1,-1,-1,-1,-1). 

Ub involves coupling between four dots, so experiences 
both quantitative three-body corrections and a four-body 
qualitative correction. Since the four spins in Ub are cou- 
pled symmetrically, the form of the four-body correction 
must also be symmetric: 



U' B 



exp 



~ E E IJ+ J' B {E AB E C D 

\A=i<j 

+EacEbd + EadEbc))] 



diag(7 3 ,7 3 ,7 3 , 

-7,-7,-7,-7, 



-7,-7,-7, 
-7,-7), 



where 7 = e~^ Js . Note that because U' B forms part of 
the Uq gate, that acts on St — states outside of the 
code space, the action of the four-body terms in it is 
nontrivial for arbitrary J' B . However, upon setting J' B to 
any integer value we recover Ub, up to an overall phase. 



4. The U 2 ,U 3 Gates 

Ui,U-i similarly involve coupling between four dots in- 
side a fixed code block, so also experience both quantita- 
tive three-body corrections and a four-body qualitative 
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correction. In this case the Heisenberg couplings are not 
symmetric, so we do not assume that the four-body terms 
are all turned on with equal coupling constants: 



Ui 



cxp 



In 



-3Eef — - (Efg + Efh + Eqh) 



+J2 E efEgh + J'2 (EacEbd + EehEfg))] 



I 5 



V P 
V 

P V* 



V 

and 

t/3 = cxp 
( 5 



, rT (- 5 + J 2+ 2 4') 

e 



-^-j= ( -3E CD - - (E AB + E AC + E BC ) 
+J' z E ab E C d + J'i (EacEbd + E ad Ebc))} 



vJ 



W P 

i7I (- 5 + J 3+ 2J 3) 

e 4^2 
P V 



where 



" = \/3(4) 2 



3(^) 2 

e = e i7r I7f( J * +2J *+ 3 ) 
( = e - i7r I75(^+ 3 ) 

/ TTV \ 

4^ COS I = 



6^4' - 164 + 1QJ k + 24 



P = 



V4V6/ 



1 , „ T ,n . / 7Tf 

) sin ■ 



(4 -4') sin 



V3l 



V4V6 

1 — e^Te 



and fc = 2 (3) for {/^ (E/3). It turns out that these two 
gates can act "classically" only when r\ = 0, which leads 



to a transcendental equation relating J' k and 4' (J'k = 
4' is one set of solutions). 



5. TTie Ui Gate 

The J7i gate is qualitatively different from all the pre- 
vious gates. It involves an interaction between five dots 
and two code blocks. In this case, the effect of the four- 
body interactions is generally to strongly interfere with 
the action of Ui . Therefore we must carefully reexamine 
this gate and consider whether it can be made compatible 
with the four-body effect. The generator of XJ\ is 



1 ° 

Hi = Ede + -z ^ Eij, 



A=i<j 

and the modified generator, in the presence of four-body 
interactions, is in general 

E 

H[=H 1 + ^ J 'ij;kl E H E kl ■ 
A=i<j<k<l 

Therefore the new gate will have the form 

There are symmetry relations between the constants 
J'ij kf th° magnitudes of the exchange interactions in H\ 
imply an equivalence between spins A, B, C, but spins 
D and E are distinct. Thus there will be four such con- 
stants, corresponding to the following sets of inequivalent 
two-body pairings: 

• one of spins {A, B, C} coupled with spin D, but 
without spin E: J' a = J' AB , CD = J' AC :BD = 

J BC;ADi 

• one of spins {A,B,C} coupled with spin E, but 
without spin D: J' b = J' AB . CE = J' AC . BE = 

J BC;AE> 

• spin D coupled with spin E: J' c = J' A b de = 
V — V 

J AC-.DE — J BC:DE' 

• one of spins {A, B, C} coupled with spin D, one 
with spin E: J d = J AD . BE = J A d-,ce = J'bd-ae 

J BD-,CE ~ J CD:AE ~ J CD;BE' 



Thus U[ can be written as 
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exp 



+EacEbd + EadEbc) 
+J'b (EabEce + EacEbe + EaeEbc) 
+J' C (EabEde + EacEde + EbcEde) 
+J'd {EadEbe + EadEce + EaeEbd 

+EaeEcd + EbdEce + EbeEcd))] 

(X+ A 

X+ A 
X+ A 

X+ A 
A X- 

A X- 
A X- 

A X- 



T_ /I 



(26) 



0/ 



where 



1/2 



V = 



x± 

A 




(9 (1 + U' h - 4J' d ) + 48 ( J' a 2 + J? -J' b J' d + J' d - J' a (A + J'd)) ) 
(3 + 8J' a 2 + 84 2 + 9f c - 4(9 + QJ' C - 4J^) + 2(3J f c - 2J' d f 
-2J' a (-3 + U' h + 3J' C - 2J' d ) - 6J' d ) 1/2 

(cos (Trs/6) ± 2iV3 {2J ' a ~ J j ~ J ' d) sin (ti/6)j (1+2^+2^+2^) 
3i(l + 2J£-2^) ^^(1+2^+2^+2^) 
e <T^(2-j;-^+2Ji) 

( cos(^) ± t(3 + 8J °- 7 j- 3J - + 2J ^ sin(^ M"^) 
V W 2V2y V V67 



/i = 



15(^-3^+2^-1) 7^ ^(2+^-2^) 

2V2iy V6 



(27) 



In comparison, Bacon's U\ gate has the form 
/ ft 



ft 



ft 



ft 



ft 



ft 



ft 



ft 



where 



ft = ie 2 V3 
H = e 73 

9 = cv. cos(-) + ^-sm(-) 



\ 
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and it may be verified that U[ reduces to U\ in the limit 
that J' a , J b , J' c , J' d -> 0. The U[ gate, like the U x gate, is 
applied at the beginning of the controlled-phase gate se- 
quence, and hence acts on the computational basis states. 
By comparing the explicit matrix representations l|26() 
and (|28|l it is clear that the crucial difference between U[ 
and U\ is the appearance of the x+ terms on the diago- 
nal of U[ (the difference between f2 and A is irrelevant: 
it translates into a global phase). In order for U[ to act 
like U\, i.e., in order for it not to prepare a superposi- 
tion of code states and the first four non-code states, x+ 
must vanish. Consulting the expression (|27fl for \+ > it is 
evident that this leads to a complicated transcendental 



equation relating the constants J' a ,J' b , J' d (but not involv- 
ing J' c ). A numerical solution of the condition x+ = 
leads to the result that the constants J' a ,J b , J' d can take 
on an infinite set of rationally related values. Upon set- 
ting the ratio J b / ' J' d to any rational number (except 1), 
there is a corresponding rational value of J' a . 



6. Summary of Conditions 

Summing up our findings, we have the following suffi- 
cient set of conditions for the modified gate sequence 



CP' = E7? 



u'MVM [u'M)u'M)\ u[(J' a , f b , J' c , J' d ) 



r 



to work as a controlled-phase gate in the presence of four- 
body interactions: 

1. The constant J' c can take on an arbitrary value. 

2. The constants J 2 and J 2 must be chosen to satisfy 
the transcendental equation r\ = 0. 

3. The constants J' 3 and J3 must be chosen to satisfy 
the transcendental equation 77 = 0. 

4. The constant J' 5 must either be zero or chosen such 

that A is an even integer (i.e. J |(Js) 2 — 2J' 5 + 1 
is an even integer). 

5. The constant J' B must be an integer. 

6. The constants J a , J b , J' d can take on an infinite set 
of rationally related values, where the ratio of any 
pair (e.g., J' b / J' d ) can be chosen completely arbi- 
trarily, and the value of the third constant is deter- 
mined by this choice. 

The most restrictive of these conditions is that J' B must 
be an integer. However, note that since the gates are ap- 
plied sequentially, this condition need only be satisfied 
during the application of the Uq gate, and it is plausible 
from the earlier sections of this paper that correspond- 
ing Heisenberg exchange constants can be found. When 
these conditions are satisfied it is indeed the case that 
CP' — (—1, 1, 1, 1) on the code space. 



E. Dimensionality of Parameter Spaces Required 
by Two-Body and Four-Body Couplings 

We caution that, although the encoding procedure de- 
scribed above has been shown mathematically to remove 



the effect of the four-body couplings, the experimental 
construction of a suitable apparatus using real quantum 
dots is another matter, as the following heuristic calcu- 
lation suggests. 

Our modified gates imply the following constraints on 
the coupling coefficients: 



U, 



(a 
(b 

U 



(d 
(e 

U. 



(f 
(g 
(h 

(i 

Ui 



(j 
(k 



gate 



K 2 [FG] =±K 2 {GH}; 

Either K 2 [FH] = 0, or k{K 2 [FH]) = 2n, where n 
must have an integer value; 

gate 



K 2 [ij] is the same for all pairs within {^4, B, C, D}; 

Ki[ABCD] = K A [ACBD] = K^ADBC]; 

K^ABCD] = 2mK 2 [AB), where m may have any 
integer value; 

gate 



K 2 [FG] = K 2 [FH] = K 2 [GH] ; 
K 2 [EF] =%K 2 [GH]; 
K 4 [EGFH] = Ki [EHFG] ; 

Either K 4 [EFGH] = K4 [EGFH] , or K 4 [EFGH] 
and K 4 [EGFH] satisfy the transcendental equa- 
tion rj{K A [EFGH], K A [EGFH]) =0; 

gate 



K 2 [AB] = K 2 [AC] = K 2 [BC}- 
K 2 [CD] =%K 2 [AB\- 
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(1) K 4 [ACBD] = Ki[ADBC]; 

(m) Either K A [ABCD] = K A [ACBD], or K A [ABCD) 
and K 4 [ACBD] satisfy the transcendental equation 
r)(K 4 [ABCD],K 4 [ACBD]) = 0; 



u{ 


gate 






(n) 


If a is the 


same for all pairs within {A, B, C, D}; 


(o) 


K 2 [DE] = 2K 2 [AB] ; 




(P) 


Ki[ABCD] 


= K 4 [ACBD] = K4ADBC]; 




(q) 


K 4 [ABCE] ■■ 


= K 4 [ACBE] = K 4 [AEBC}; 




w 


K A [ADBE] 


= K 4 [AEBD] ; 




(s) 


Ki[ADCE] 


= Ki[AECD\; 




ft) 


Ki[BDCE] 


= K 4 [BECD}; 




(u) 


Ki[ADBE] 


= Ki[ADCE] = K 4 [BDCE\- 




(v) 


Ki[ABDE\ 
K 4 [ADBE\; 


is a single-valued function 


of 


(w) 


Ki[BCDE] 
Ki[ADCE]\ 


is a single-valued function 


of 


(x) 


Ki[ACDE] 
K 4 [BDCE}. 


is a single-valued function 


of 



Since the coupling coefficients must in general vary 
with time in order to satisfy all of these constraints (for 
example, K 2 [CD] and K 2 [AC] would be equal during the 
operation of U' B , but unequal during U^), we also assume 
that particular constraints need to be concurrently satis- 
fied only when they arise from the same gate. 

First, by the same reasoning used to derive Eqs. © 
and l|18|l . we note that a four-dot Hamiltonian for the 
geometry of {A, B, C, E} contains a constant term and 9 
independent coupling coefficients. If these 9 coefficients 
take on a given set of values and we wish to adjust them 
to meet constraints such as those listed above, we would 
need 9 additional degrees of freedom in the system. We 
make the conservative assumption, however, that one 
two-body coefficient and one four-body coefficient can 
be left unaltered and the others adjusted to correspond 
to them, which means that only 7 additional parame- 
ters are required. Similarly, for the subset {A, B, D, E} 
({A,C,D,E}, {B,C,D,E}, {A,B,C,D}), there are 9 
(6, 7, 5) independent coupling coefficients, for which we 
require 7 (4, 5, 3) tunable parameters if a given set of con- 
straints are to be satisfied. Of course, we will count one 
more degree of freedom whenever a constraint includes 
relationships between the two-body and four-body ener- 
gies. 

Now suppose that we designate one "base" choice of 
{xb, x c , x v }, such that within each of the two squares, all 
the quantities K 2 [ij] are equal, all the quantities K 4 [ijkl] 
are equal, and K 4 [ijkl] — 2K 2 [ij]. That arrangement 



can simultaneously satisfy constraints (b), (c), (d), (e), 

( f )> ( h ), (i)> 0), Q)> ( m ). ( n )> and (P). provided that 
the value of K 2 [FH] is chosen appropriately. From this 
potential, we would need to make one change within 
{E, F,G, H} to reach condition (a) or condition (g), or 
one change within {5,6,7,8} to obtain (k) or (o). The 
couplings of {A, B, C, E} must be adjusted to match (q) 
while still satisfying (n), (o), and (p), which requires 6 
additional degrees of freedom as explained in the previ- 
ous paragraph. Similarly, (v) ((w), (x)) and (r) ((s), (t)) 
together imply particular adjustments to the four-body 
couplings in {A, B, D, E} ({A,C,D,E}, {B,C,D,E}), 
which requires 6 (3,4) new parameters. (The single- 
valued function in question is the same for all three 
cases, so K 4 [ABDE] ends up equalling K 4 [ACDE] and 
K 4 [BCDE\.) Finally, we need two more degrees of free- 
dom available somewhere in order to meet constraint (u) , 
for a grand total of 28 degrees of freedom. 

To put the size of this number into perspective, we will 
also count the independently tuned energies necessary 
to meet the conditions on EQC using pairwise couplings 
alone. By choosing a suitable combination of {xb, x c , x v } 
for an entire eight-spin system, we could satisfy (b), (c), 
(f), (j), and (n) at the same time; one more degree of free- 
dom would be necessary to also satisfy (o). Starting from 
such a system, we could presumably satisfy (a) or (g) 
by adjusting one parameter within {E, F,G, if}, or sat- 
isfy (k) by adjusting one parameter within {A, B, C, D}. 
Hence we estimate that 7 degrees of freedom are required 
for the purely Heisenberg Hamiltonian used in Ref. [30| . 
We see that, even if one presupposes the ability to create 
and position many identical qubits of the form 11' I'D , i|23|) 
(3 free parameters), accounting correctly for two-body 
and four-body coupling is still a great deal more demand- 
ing than two-body coupling alone. It is this experimen- 
tal challenge that must be weighed against the increased 
length (and hence vulnerability to decoherence) of pulse 
sequences employing only two-body couplings |3l| . 



V. SUMMARY AND CONCLUSIONS 

Earlier work [l7l IT^ | showed that in highly symmet- 
rical geometries, the interaction between three and four 
mutually interacting electrons confined in parabolic po- 
tentials contains many-body terms, which in the case of 
four electrons qualitatively modify the usual Heisenberg 
interaction. In this work we have improved upon these 
early results by considering realistic, linear and square 
geometries, and by utilizing Gaussian confining poten- 
tials. Specifically, we have shown in a Heitler-London 
calculation that in the case of four mutually interacting 
electrons, in both the linear and square geometries, the 
system's Hamiltonian contains four-body exchange terms 
which may be of comparable strength to the Heisenberg 
exchange interactions. This can have important impli- 
cations for quantum information processing using cou- 
pled quantum dots. We have considered, in particular, 
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the implications for the quantum computing using log- 
ical qubits encoded into decoherence-free subspaces of 
four electrons per qubit. We showed that previously de- 
signed conditional quantum logic gates between these en- 
coded qubits must be modified, in order to account for 
the four-body terms that alter the (previously assumed) 
Heisenberg interaction, when four or more electrons are 
coupled simultaneously. This requires the ability to tune, 
to a certain extent, the four-body exchange constants. It 
is worth noting, however, that there are alternatives to 
this method of implementing encoded conditional logic 
gates, which may be less demanding. In particular, it is 
worth exploring the possibility of completing the set of 
universal encoded quantum logic gates by supplementing 
single-qubit gates (where, as we have shown, four-body 



effects are harmless) with measurements and teleporta- 
tion, as in linear optics quantum computing |4f| . This 
will be a subject for future research. 
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FIG. 1: (Color online) Plot of Ko, the overall energy shift, as 
a function of dimensionless barrier height X), and overall well 
depth x v in the case of three mutually interacting electrons in 
a linear geometry. In this and succeeding figures, the Coulomb 
repulsion parameter x c is set to 1.5 as in Ref. |17j . 




FIG. 2: (Color online) Plot of K 2 [AB], the two-body coupling 
coefficient for adjacent dots, as a function of dimensionless 
barrier height xj, and overall well depth x v in the case of 
three mutually interacting electrons in a linear geometry. 



FIG. 3: (Color online) Plot of K2[AC], the two-body coupling 
coefficient for non-adjacent dots, as a function of dimension- 
less barrier height Xb and overall well depth x v in the case of 
three mutually interacting electrons in a linear geometry. 




FIG. 4: (Color online) Plot of Ko, the overall energy shift, as 
a function of dimensionless barrier height Xb and overall well 
depth x v in the case of four mutually interacting electrons in 
a square geometry. 



FIG. 5: (Color online) Plot of K2[AB], the two-body coupling 
coefficient for adjacent dots, as a function of dimensionless 
barrier height Xb and overall well depth x v in the case of four 
mutually interacting electrons in a square geometry. 
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FIG. 6: (Color online) Plot of K 2 [AC], the two-body coupling 
coefficient for non-adjacent dots, as a function of dimension- 
less barrier height x\, and overall well depth x v in the case of 
four mutually interacting electrons in a square geometry. 



FIG. 7: (Color online) Plot of K 4 [ABCD], the four-body 
coupling coefficient for pairs of adjacent dots, as a function 
of dimensionless barrier height Xb and overall well depth x v 
in the case of four mutually interacting electrons in a square 
geometry. 




FIG. 8: (Color online) Plot of K 4 [ACBD], the four-body cou- 
pling coefficient for pairs of non-adjacent dots, as a function 
of dimensionless barrier height Xb and overall well depth x v 
in the case of four mutually interacting electrons in a square 
geometry. Note that two of the axis directions are reversed 
from the preceding figures. 
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FIG. 9: (Color online) Partitioning of the Hilbert space of 
N spin- 1/2 particles into DF subspaces (nodes of the graph). 
The integer above each node represents the number of paths 
leading from the origin to that node. 



